Topological superfluid in one-dimensional spin-orbit coupled atomic Fermi gases 
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We investigate theoretically the prospect of realizing a topological superfluid in one- dimensional 
spin-orbit coupled atomic Fermi gases under Zeeman field in harmonic traps. In the absence of 
spin-orbit coupling, it is well-known that the system is either a Bardeen-Cooper-Schrieffer (BCS) 
superfluid or an inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superfluid. Here we show 
that with spin-orbit coupling it could be driven into a topological superfluid, which supports zero- 
energy Majorana modes. However, in the weakly interacting regime the spin-orbit coupling does 
not favor the spatially oscillating FFLO order parameter. As a result, it seems difficult to create an 
inhomogeneous topological superfluid in current cold-atom experiments. 
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I. INTRODUCTION 

Topological superfluids are new states of matter that 
attract intense attentions in recent years P, 01 ■ They 
have a full pairing gap in the bulk and exotic gapless ex- 
citations at the edge - the so-called Majorana fermions 
- which obey non-Abelian statistics [3, 4]. These excita- 
tions are immune to decoherence caused by local per- 
turbations. By properly braiding excitation quasipar- 
ticles, topological quantum information might be pro- 
cessed. As a result, topological superfluids could provide 
an ideal platform for topological quantum computation 

Q . Because of this potential application, the realiza- 
tion of topological superfluids in a well-controlled envi- 
ronment is highly desirable. 

Theoretically, there are a number of proposals on re- 
alizing a topological superfluid in two-dimensional (2D) 
settings, including the use of 2D p-wave pairing 0, Q, 
proximity coupling to a conventional s-wave supercon- 
ductors for the surface state of three-dimensional (3D) 
topological insulators [E^JT], and 2D atomic Fermi gases 
with strong Rashba spin-orbit coupling (T2| . ITsj . It is 
also possible to create a topological superfluid in one- 
dimensional (ID) solid-state systems by suitably engi- 
neering spin-orbit coupling of electrons, such as InAs 
wires and banded carbon nanotubes |14l - [l7| . The pur- 
pose of this work is to examine the possibility of observ- 
ing topological superfluids in ID ultracold atomic Fermi 
gases lis], which may be regarded as highly controllable 
quantum simulators of the corresponding ID solid-state 
systems. We note that ID atomic Fermi gases can now 
be routinely created in cold-atom laboratories [19]. The 
spin-orbit coupling for neutral atoms may also be gener- 
ated by using the so-called "non-Abelian synthetic gauge 
fields" technique [20, 21]. 

Even in the absence of spin-orbit coupling the ID ul- 
tracold atomic Fermi gas is of great interest. It hosts 
a Bardeen-Cooper-Schrieffer (BCS) superfluid and an 
exotic inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) superfluid [H, [H-ii], respectively, in the case 



of balanced and imbalanced spin-populations. Here we 
show that by adding spin-orbit coupling both superflu- 
ids can turn into a topological superfluid. We discuss 
in detail the resulting zero-energy Majorana edge modes 
and their possible experimental signature. We also ex- 
plore the possibility of creating an inhomogeneous topo- 
logical superfluid with spatially oscillating FFLO order 
parameter. Unfortunately, the spin-orbit coupling seems 
to suppress the FFLO order parameter. As a result, in 
the weakly interacting regime we always find the same 
topological superfluid with a uniform order parameter, 
whatever the initial state is a BCS or FFLO super- 
fluid. Our study is based on the self-consistent solution of 
fully microscopic Bogoliubov-de Gennes (BdG) equations 
[25I . [28| . It enables ab-initio simulations under realistic 
experimental conditions. 

The paper is organized as follows. In the next section 
(Sec. II), we present the model Hamiltonian and the BdG 
equations. In Sec. Ill, we discuss the phase diagram at a 
sufficiently large spin-orbit coupling and the phase tran- 
sition from BCS superfluid to topological superfluid. The 
wave-functions of Majorana edge modes are shown and 
their possible experimental detection is considered. In 
Sec. IV, we present the phase diagram at a given Zee- 
man fleld and show the transition from FFLO superfluid 
to topological superfluid. Finally, in Sec. V we provide 
conclusions and some flnal remarks. 



II. MODEL HAMILTONIAN AND BDG 
EQUATIONS 

We consider a trapped two-component ID atomic 
Fermi gas under a non-Abelian gauge fleld (spin-orbit 
coupling) and Zeeman field, described by the model 
Hamiltonian, 



H 



J dxip'^ (x) [Hq (x) - hoz + Xkuy] (x) 
+91D / dxi)l{x)i)l{x)i)i{x)i}^{x) , (1) 
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where tp'^ (x) = (x) , -01 (x)] denotes collectively the 
creation field operators for spin-up and spin-down atoms. 
In the single-particle Hamiltonian (i.e.. the first line 
of the above equation), "Hg (x) = —{h^/2m)d'^/dx'^ + 
raiiP'x^ 1^ — /i describes the single-particle motion in a 
harmonic trapping potential muj^x'^ /2 and in reference to 
the chemical potential /j,, the strength of the Zeeman field 
is denoted by /i, Xkuy = ~iX{d/dx)ay is the spin-orbit 
coupling term with coupling strength A, ay and ct^ are 
the 2x2 Pauli matrices. The second line of the equation 
is the interaction Hamiltonian, with the (attractive) in- 
teraction strength given by the s-wave scattering length: 
giD = -2h^ /{maiD). 

The model Hamiltonian Eq. ([T]) can be realized 
straightforwardly with cold fermionic atoms. It is a di- 
rect generalization of the standard model Hamiltonian 
for a ID spin-imbalanced Fermi gas, through the inclu- 
sion of a non- Abelian synthetic gauge field Xkuy . Exper- 
imentally, a bundle of ID spin-imbalanced atomic Fermi 



gases can now be manipulated using 2D optical lattices 
|19| . The generalization of the synthetic gauge field Xkuy 
has already been demonstrated in a 3D Bose gas of *^Rb 
atoms fSo'l . In addition, its realization in fermionic atoms 
has been proposed [21]. Therefore, all the techniques re- 
quired to simulate Eq. (1) are within current experimen- 
tal reach. 

To understand the ID superfiuidity in the presence of 
spin-orbit coupling, we calculate elementary excitations 
within the mean- field BdG approach jMUl- The wave- 
function of low-energy fermionic quasiparticles VP,, {x) 
with energy i?^ is solved by, 

•HsdG*^ {x) = E^n'ir, {x) , (2) 

where {x) = [u^y^ (x) , u^r, {x) , wt») {x) , vir, {x)]'^ in the 
Nambu spinor representation and the BdG Hamiltonian 
'Hsdc reads accordingly. 



T-LbcIG 



n^{x) - h -Xd/dx 
Xd/dx n^{x) + h 

A*(a;) 
-A*(x) 



-A(a;) 
A(.t) 
-n^{x) + h Xd/dx 
~Xd/dx -ni{x) - h 



(3) 



Here A(x) = -(5id/2) E,K.^%/(-E^r,) + 
uirii>^^f{—Eri)] is the order parameter and 
f {x) = l/[e'^/('^^-^^+] is the Fermi distribution 
function at temperature T. The order parame- 
ter is to be solved self-consistently together with 
the number equation for the chemical potential, 
J drln-f (r) -t- (r)] = N, where TV is the total number 
of atoms and the density of spin-cr atoms is given 
by, n,ix) - il/2)j:,[Kn\' fiErj) + Krf f{~E,)]. 
We note that the use of Nambu spinor representation 
leads to an inherent redundancy built into the BdG 
Hamiltonian 0. HbcIG is invariant under the particle- 
hole transformation: (x) — >■ v* (x) and E^i — >■ —E^. 
Therefore, every eigenstate with energy E has a partner 
at —E. These two states describe the same physical 
degrees of freedom, as the Bogoliubov quasiparticle 
operators associated with them satisfy F^ = fI^. This 
redundancy has been removed by multiplying a factor of 
1/2 in the expressions for order parameter and atomic 
density. 

The BdG equation ^ can be solved by expanding 
Uari (x) and Vari (x) in the basis of ID harmonic oscil- 
lators. On such a basis, Eq. ^ is converted to a secu- 
lar matrix. A matrix diagonalization then gives the de- 
sired quasiparticle energy spectrum and wave-functions. 
Numerically, we have to truncate the summation over 
the energy levels rj. For this purpose, we adopt a hy- 
brid strategy developed earlier by us for an imbalanced 



Fermi gas without spin-orbit coupling [2^ [2^. We in- 
troduce a high energy cut-off Ec, above which a lo- 
cal density approximation (LDA) is used for the high- 
lying energies and wave-functions. This leads to an ef- 
fective coupling constant in the gap equation, A (a;) = 

-[gig (x) /2] Eni^tvvUfiEv) + Mi,i'|,/(-i?,)], where 
J2ri is now restricted to | < Ec- We refer to Ref. 

f^Bj for further details of (x) and the LDA atomic 
density. 

In harmonic traps, it is useful to characterize the in- 
teraction strength by using a dimensionless interaction 
parameter [2^, 7 = —mgi]j/{h'^no) = 2/(noai£)), where 
riQ is the zero-temperature center density of an ideal 
two-component Fermi gas with equal spin populations 
A^/2. In the Thomas-Fermi approximation (or LDA), 
no = 2N^/'^ /{iraho) and a^o = ^Jh/{mui) is the char- 
acteristic oscillator length of the trap. Therefore, the 
dimensionless interaction parameter is given by, 

T = (— ) ■ (4) 

We note that, for a ID atomic Fermi gas created us- 
ing 2D optical lattices, the typical dimensionless interac- 
tion strength is about 7 = 3 ~ 5 [l9l [25| . Through- 
out the paper, we shall take a slightly smaller value 
of 7 = 7r/2 ~ 1.6, in order to validate the mean-field 
treatment. It is also convenient to use the Thomas- 
Fermi energy Ep — {N/2)huj and Thomas-Fermi radius 
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xp = N^^'^O'ho as the units for energy and length, re- 
spectively. For the spin-orbit coupling, we use a dimen- 
sionless parameter Xkp/Ep, where kp = ^/2mEp is the 
Thomas- Fermi wavevector. We have performed numeri- 
cal calculations for a Fermi gas of = 100 fermions in 
traps at both zero temperature and finite temperature. 
In the following, we present only the zero-temperature 
results, as the inclusion of a finite but small temperature 
(i.e., T — O.lTp) essentially does not affect the results. 
The Fermi energy is Ep = {N/2)huj = 50huj. We have 
taken a cut-off energy Ec = 'iEp = 200hijj and have used 
37V = 300 ID harmonic oscillators as the expansion func- 
tions. These parameters are already sufficiently large to 
ensure the accuracy of calculations. 

III. PHASE DIAGRAM AT A GIVEN 
SPIN-ORBIT COUPLING 

The most salient feature of a spin-orbit coupled Fermi 
gas is the appearance of topological superfluidity and 
zero-energy Majorana fermion mode, under an appropri- 
ate Zeeman field. The quasiparticle operators of Ma- 
jorana fermions are real and satisfy 7 = 7^, which 
means that a quasiparticle is its own antiparticle 0, 0| . 
Mathematically, we can always write a complex ordi- 
nary fermion operator c in terms of two real Majorana 
fermions 71 and 72, such as c = 71 — 172. An ordinary 
fermion may therefore be viewed as a bound state of 
two Majorana fermions, which in general can not be de- 
confined. However, the deconfinement does happen in a 
topological superfluid, leading to two Majorana fermions 
localized respectively at the two edges of topological su- 
perfluid. This can be clearly seen with the help of the 
particle-hole redundancy of the BdG equation 0, hj. Let 
us image that we have a zero-energy solution E = 0. Be- 
cause of the particle- hole redundancy F^; = F^^;, we will 

immediately have Fq = fJ - exactly the defining feature 
of a Majorana fermion. We note that, zero-energy Ma- 
jorana fermions should always come in pairs, since the 
original model Hamiltonian describes ordinary fermions 
only and each Majorana fermion is just a half of ordi- 
nary fermion. It is straightforward to check from the 
BdG Hamiltonian that the wave functions of two paired 
Majorana fermions should satisfy Ua- (x) = v* (x) and 
Ua (x) = —v*{x), respectively. The former follows the 
particle-hole symmetry, while the later is required to ex- 
press an ordinary fermion by two Majorana fermions (3^ . 

To satisfy the prescription of a zero-energy solution for 
a topological superfiuid, the quasiparticle energy spec- 
trum must become gapless at a certain point. In the 
case of a homogeneous spin-orbit coupled Fermi gas un- 
der a Zeeman field, this happens at a critical Zeeman 
field 

The system will be in a conventional superfluid at h < he 
and in a topological superfluid aX h > he- For a trapped 



system, however, the critical Zeeman field may become 
position dependent. As a result, in harmonic traps we 
would have a mixed phase with both conventional and 
topological superfluid components, which separate spa- 
tially in real space. Without confusion, we shall still 
refer to such a mixed phase as a topological superfluid. 

A. Phase diagram at Xkr/Ep = 1 
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Figure 1: (color online) Phase diagram at a given spin-orbit 
coupling Xkr/Ep = 1, determined from the behavior of 
the lowest eigenenergy of Bogohubov quasiparticle spectrum, 
min{|£',,|}. As the Zeeman field increases, the system evolves 
from a conventional BCS superfluid to a topological super- 
fluid, and finally to a normal state. The two insets in the mid- 
dle and right show the quasiparticle spectrum at h/ Ef = 0.3 
and 0.5, respectively. 

In Fig. 1 we report the phase diagram at a fixed spin- 
orbit coupling strength Xkp/Ep — 1. The emergence of 
a topological superfluid can be clearly revealed by the 
behavior of the lowest eigenenergy of the quasiparticle 
energy spectrum. As shown in the middle inset, at a 
small Zeeman field the energy spectrum is gapped. How- 
ever, by increasing the Zeeman field above a critical value 
of h ^ 0.35Ep, the lowest eigenenergy becomes exponen- 
tially small. Four quasiparticle modes with nearly zero 
energy appear, as seen clearly from the right inset. By 
further increasing the Zeeman field {h > 0.65Ep), the 
system will be driven into a normal state with negligible 
superfiuid order parameter. 

The appearance of the topological superfluid can also 
be monitored by the calculation of h — hc{x), where 
hc{x) = ^ /i^(x) -1- A2(x) is the local critical Zeeman held 
for a local uniform cell at position x with the local chem- 
ical potential /i(x) = ^ — muP'x^ jl and order parameter 
A(a;). The local uniform cell would be in the topologi- 
cal superfluid state if ft. > hdx). In Fig. 2, we present 
h — hc{x) and A(a;) at different phases. In accord with 
Fig. 1, at a smaU field h = 0.3Ep (Fig. 2(a)), h < hc{x) 
for any position x and the whole Fermi cloud is in the con- 
ventional superfluid. At the held h = O.bEp (Fig. 2(b)), 
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Figure 2: (color online) Spatial dependence of the critical 
Zeeman field h — hc{x) (solid lines) and the superfluid order 
parameter A(x) (dot-dashed lines), at Xkr/Ep = 1 and at 
three Zeeman fields h/Ep ~ 0.3, 0.5, and 0.8. The cross- 
patterns highlight the area in which the atoms are in the 
topological superfiuid state. 



we find h > hc{x) at the two wings of the harmonic trap 
and therefore there are two blocks of topological super- 
fluid, as highlighted by the cross-pattern. At an even 
large Zeenian-ficld (Fig. 2(c)), the area of h > hc{x) 
extends over the whole system. However, the superfluid 
order parameter becomes so small, the system can no 
longer be viewed a superfluid. We note that, at large 
attractive interactions where the order parameter is not 
destroyed by large Zeeman field, it is possible to have a 
single topological superfluid throughout the whole Fermi 
cloud. 



B. Majorana fermions 

In each of the topological superfluid phases, we should 
flnd two Majorana fermion modes, well-localized at the 
two edges respectively. At the Zeeman fleld h = 0.5Ep, 
we therefore could have four Majorana fermions, as in- 
dicated by the energy spectrum in the right inset of Fig. 
1. The wave functions of these Majorana fermions are 
shown in Figs. 3 and 4 for states localized at a; ~ ±0.5xi? 
and ±l.la;i?, respectively. It is interesting that the wave 
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Figure 3: (color online) Wave functions of the paired Majo- 
rana modes at the inner wing of the trap, x ~ zbO.Sa;^: one 
has the energy Ezes — —3.3 x 10~^_Ef (a) and the other 
EzES — +3.3 X Ef (b). Both modes satisfy the symme- 
try requirement for Majorana wave-functions. The wave func- 
tions are in units of a'^y^ ■ Here h = O.S-Ef and Xkr/Ep ~ 1. 
The system is in the topological superfluid phase. 




Figure 4: (color online) Wave functions of the paired Majo- 
rana modes at the outer wing of the trap, x ~ ±1.13;^: one 
has the energy Ezes — —7.2 x 10~^''iJF (a) and the other 



Ezes ^ +7.2 x 10" 
as in Fig. 3. 



Ef (b). Other parameters are the same 
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functions of two paired Majorana fermions, for example, 
these located at x ~ —O.bxp and x ~ +0.5a;F (Fig. 3), 
tend to interfere with each other [l^, [l^- This quasi- 
particle interference or tunneling leads to the splitting of 
degenerate zero energy Majorana modes to a finite but 
exponentially small energy: Ezes — ±3.3 x IQ^^Ep- 
The tunneling between the paired Majorana fermions at 
the outer wing of the trap, x ~ zbl.lxF, is more diffi- 
cult (see Fig. 4), so the energy splitting is much smaller, 
i.e., Ezes — ±7.2 x lQ~^'^Ep. It is readily seen that the 
paired wave functions satisfy either u„ [x) — v* (x) or 
Ua {x) = —V* (x), as anticipated by the required symme- 
try of Majorana wave functions. 



C. Density distribution and local density of states 

We now consider the possible experimental signature 
for observing topological superfluid and the associated 
Majorana fermions. The useful experimental tools in- 
clude in-situ absorption imaging and spatially resolved 
radio- frequency (rf) spectroscopy [30], which give respec- 
tively the density distribution and the local density of 
states of the Fermi cloud [3l|. 

In Fig. 5, we plot the spin-up n-f{x) and spin-down 
714.(2;) density distribution and their difference An{x) = 
'nf{x) — n,i{x) at different phases. While the shape 
of the spin- up density distribution n^{x) is nearly un- 
changed across different phases, in the topological super- 
fluid phase (see Fig. 5(b) at h = 0.5Ep) the spin-down 
density distribution n^ix) shows an interesting bi-modal 
structure. It decreases rapidly when the atoms enter the 
topological area from the center. Accordingly, a broad 
dip appear in the density difference around the trap cen- 
ter. The bi-modal distribution in n^^x) may be regarded 
as a useful and convenient feature to identify the topolog- 
ical superfluid. However, it is not a characteristic feature 
for identifying the Majorana modes, as the contribution 
of the Majorana modes to the density distribution is neg- 
ligibly small, i.e., relatively at the order of N~'^/^. 

A practical way to probe the Majorana fermions is to 
measure the local density of states using the spatially 
resolved rf spectroscopy [30', 'ST], with which we antici- 
pate that the contributions of Majorana fermions will be 
well-isolated in both energy domain and real space. The 
local density of states for spin-up and spin-down atoms 
is defined by. 



Pa{x,E) = ^Y1 

V 



^arj I 



■5{E~E,-,) + 



■5{E + E,,) 



(6) 

In Fig. 6, we report the local density of states in the topo- 
logical superfluid state. Near the zero energy, the contri- 
butions from Majorana fermions are clearly visible and 
are well-separated from other quasiparticle contributions 
by an energy gap A ~ Q.lEp. It is interesting to note 
that the Majorana modes at a; ~ ±1.10:^ and ±0.5xf 
contribute to p-^{x,E) and pi{x^E), respectively. This 
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Figure 5: (color online) The spin-up and spin-down density 
distribution, n-|-(a;) (dashed lines) and n^(a;) (solid lines), and 
their difference An(a;) = nt(^) ~ '"■li?^') (dot-dashed lines), 
are shown at the conventional superfluid phase (a), topolog- 
ical superfluid state (b), and normal state (c). The den- 
sity distributions are in units of the Thomas-Fermi density 
no = 2N^^^ / ijiaho)- The spin-orbit coupling is Xkr/Ep = 1. 




Figure 6: (color online) Linear contour plot (in arbitrary 
units) for the local density of states of spin-up atoms p-f{x, E) 
(a) and of spin-down atoms pi(x,E) (b). The contributions 
from Majorana fermions are highlighted by circles. Here the 
Fermi cloud is in the topological superfluid state with pa- 
rameters h = 0.5Ef and Xkp/EF ~ 1. In the calculation, 
the (5-function in pa{x,E) has been simulated by a Lorentz 
distribution with a small energy broadening F = O.OlEp. 
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can be understood from the wave function of Majorana 
modes, as shown in Figs. 3 and 4. The wave-functions 
at ±0.5xf are dominated by the spin-down component, 
while the wave- functions at ±1.1xf have mainly the spin- 
up component. 



IV. PHASE DIAGRAM AT A GIVEN ZEEMAN 
FIELD 

We now turn to consider the possibility of observing 
a topological superfluid with spatially oscillating order 
parameter [23, [2^. In the absence of spin-orbit coupling, 
it is known that the ground state of an imbalance ID 
Fermi gas under Zeeman field can be an inhomogeneous 
FFLO superfluid with oscillating order parameter. It is 
therefore natural to ask: what is the fate of such a FFLO 
superfluid when we switch on the spin-orbit coupling? 
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Figure 7: (color online) Phase diagram at a given Zeeman 
field h — OAEf, determined from the behavior of the lowest 
eigenenergy of Bogoliubov quasiparticle spectrum, min{|iJ^|}. 
As the spin-orbit coupling increases, the system evolves from a 
FFLO superfluid to a topological superfluid. The inset shows 
the critical Zeeman field h — hc{x) and the order parameter 
A(a;) at Xkr/Ep — 0.3, where the Fermi gas is in the FFLO 
superfiuid state. 

In Fig. 7, we present the phase diagram at a given 
Zeeman field h = OAEp, determined again by tracing 
the behavior of the lowest eigenenergy of the quasiparti- 
cle spectrum as a function of the spin-orbit coupling. The 
density distributions and order parameter are reported in 
Fig. 8 for three values of spin-orbit coupling. At small 
spin-orbit coupling, we find a stable FFLO order parame- 
ter which is modified slightly by the spin-orbit coupling. 
However, in the area where App^o (x) is nonzero, the 
criterion for a topological superfluid h > hc{x) is al- 
ways not satisfied, as seen from the inset for the case of 
Xkp/Ep = 0.3. This excludes the coexistence of FFLO 
superfluid and topological order. As a result, the energy 
spectrum is gapped and min{|i?,J} > 0. With increas- 
ing the spin-orbit coupling above Xkp/Ep ~ 0.6, we ob- 
serve that the lowest eigenenergy becomes exponentially 



0.6 



< 0.4 



i 0.2 



0.0 



(a) U/£^ = 0.2 




(b) XkJE^ = 0.6 

















Figure 8: (color online) Density distributions and order pa- 
rameter at /i = QAEf and at three different spin-orbit cou- 
plings: Xkp/EF = 0.2 (a), 0.6 (b), and f.O (c). The density 
distribution Uc^ix) is in units of the Thomas- Fermi density 
no = 2N^^'^ /{naho)- The order parameter A(a;) is in units of 
the Fermi energy Ep. 



small, suggesting a topological superfluid. However, in 
this case, the order parameter no longer oscillates in real 
space, as shown in Figs. 8(b) and 8(c). Therefore, we 
conclude that it seems impossible to create an inhomo- 
geneous topological superfluid with spatially oscillating 
order parameter in ID spin-orbit coupled Fermi gas, if 
we do not tailor specifically the geometry or other pa- 
rameters of the Fermi cloud. 



CONCLUSIONS 



In conclusions, we have investigated theoretically the 
properties of a ID imbalanced Fermi gas under non- 
Abelian synthetic gauge field. We have predicted that 
by suitably tuning the strength of spin-orbit coupling and 
Zeeman field, it is possible to create a topological super- 
fluid, which hosts Majorana zero-energy fermions at its 
edge. The order parameter in the topological superfluid 
is always of the conventional Bardeen-Cooper-Schrieffer 
type, as the spin-orbit coupling tends to destroy inhomo- 
geneous Fulde-Ferrell-Larkin-Ovchinnikov pairing. Ex- 
perimentally, the topological superfluid may be identified 
from the bimodal distribution of the spin-down atomic 
density by using in-situ absorption imaging. The as- 
sociated Majorana fermions may be detected by apply- 
ing the spatially resolved radio-frequency spectroscopy, 
which would show a well-isolated signal at zero energy. 

At the end of this paper, we would like to emphasize 
that the ultracold atomic Fermi gas with non-Abelian 
synthetic gauge field is an ideal platform for creat- 
ing topological superfluid and manipulating Majorana 
fermions, because of its unprecedented controllability and 
flexibility. This system can now be readily realized in ul- 
tracold atom laboratories. 
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